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Abstract: We construct the classical and canonically quantized theories of a massless 
scalar field on a background lattice in which the number of points — and hence the number 
of modes — may grow in time. To obtain a well-defined theory certain restrictions must be 
imposed on the lattice. Growth-induced particle creation is studied in a two-dimensional 
example. The results suggest that local mode birth of this sort injects too much energy 
into the vacuum to be a viable model of cosmological mode birth. 



Keywords: sos, mqg, lqf, lmg. 



Contents 



[l]. Introduction 

|l.l| Overview of the approach and results 



I 



|2]. Latt ice structure 

.1| Hyper diamond lattice 



I 



Lattice dynamics 

Hyperdiamond lattice action 



3.0.1 



3.1 



3.2 



3.4 



3.0.2 Irregular lattice action 
Discrete field equations 
Initial value problem 
Phase space of half-solutions 
Symplectic form 

Degeneracy 



3.4.1 



3.5 



Poisson brackets 

Local fields and locality 



3.5.1 



Quantization 

QFT on a hyperdiamond lattice 

Hyperdiamond modes 
Quantization on a hyperdiamond lattice 
Harmonic modes 
Zero-modes 
5.2.31 Field operator and Fock vacuum 



5.2.1 



5.2.2 



Mode birth on a growing two-dimensional lattice 
|6.1| Regular diamond lattice 



6.2 Growing lattice: mode birth 



6.3 Mode decompositions and Fock space 

6.3.1 Evolution rule 

6.3.2 Choice of birth modes 



6.4 Growth-induced particle creation 



nscussion 



S Dis 
\A[ Bogoliubov coefficients 








10 
10 

n 

U 

n 

w 

16 
17 
17 

18 

19 
19 

m 
m 

M 

24 
25 



- 1 - 



1. Introduction 



If, as suspected by many, there is a short distance cutoff on physical degrees of freedom, then 
as the universe expands the number of degrees of freedom must be growing. Such a scenario 
might neatly account for the otherwise puzzlingly low initial entropy of the universe, and 
at long wavelengths might be indistinguishable from the usual continuum physics in an 
expanding background spacetime. But to explore this avenue is difficult, since it raises the 
problem of formulating dynamics in a setting where the number of degrees of freedom is 
not the same at each time [[ij. 

In the framework of effective field theory (EFT), it is straightforward to describe mode 
birth. Consider by analogy a growing droplet of liquid. The collective vibrational modes 
define an EFT with a cutoff at the atomic scale. As new atoms condense onto the droplet, 
new modes are added; but at the growing scale, the EFT does not apply. One can work 
in the Fourier transform domain, not with a distance cutoff but with a wavevector cutoff. 
Then EFT modes can be inserted by hand in the adiabatic ground state, but in this 
approach one learns nothing about the microscopic physics associated with growth of the 
microstructure of spacetime. 

Perhaps a quantum gravity theory can be formulated in which states with different 
numbers of degrees of freedom are all contained in the same Hilbert space, so that in the pre- 
geometric sense the number of degrees of freedom is actually constant. Then presumably 
mode birth would be just a feature of the dynamics. This is a deep problem. Here we 
try to gain some insight by considering what should be a much simpler problem: quantum 
field theory (QFT) on a background growing, or more generally irregular lattice. 

An appealing idea for a fundamental discrete, irregular structure is that of a causal 
set 0, §]. It is an outstanding problem whether and how one might formulate quantum 
field theory on a causal set. (Some ideas on classical field theory have been worked out 
however Q.) This question forms a second motivation for the present work, although the 
causal sets we consider have two properties not shared by a generic causal set: finite linking 
and a 'layered' structure. 

We thus wish to consider a simple free scalar field theory model on an irregular lattice 
structure. The challenge is to formulate a theory that allows for the increase of degrees of 
freedom, while possessing a good quantum mechanical interpretation and agreeing at long 
wavelengths with standard quantum field theory in an expanding background spacetime. 

1.1 Overview of the approach and results 

It is by no means clear how to modify quantum field theory to allow for irregular lattices, 
in particular growing ones. One could imagine just writing down an action, and a cor- 
responding path integral. This has been done on random lattices in Euclidean signature, 
and the path integral can be given there a statistical probability interpretation. But in 
the Lorentzian context, without the usual underlying Hamiltonian structure there is no 
guarantee that the complex amplitudes so defined would have any consistent probability 
interpretation at all. For example the transition amplitudes would in no sense be unitary. 
One could try to use the decoherence functional approach, rather than demanding strict 
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unitarity. In fact just such an approach was explored in || in a very general setting of evo- 
lution on causal networks. But there is no guarantee that useful sets of decohering histories 
could be found in such a general setting. Abstract formulations of quantum dynamics on 
causal sets have also been developed || but without direct relation to the dynamics of 
any particular field theory. For the present work the aim is simply to recover the quantum 
field theory of a linear field in the large, so we prefer to stay closer to a local Hamiltonian 
formulation, adopting the algebraic approach to QFT||, M. It may be that for the class 
of lattices considered here, the path integral approach does make sense, but we have not 
investigated that question. 

The algebraic approach has for ingredients a C*-algebra and a state. For example, 
the algebra might be generated by the abstract quantum field "operator" satisfying a field 
equation and canonical commutation relations, and the state is a positive linear functional 
on the algebra that assigns expectation values to observables. The conservation of a sym- 
plectic structure under evolution from one time slice to the next plays a critical role in the 
way both of these structures are usually defined in practice. The equation of motion must 
be compatible with the commutation relations, a property whose validity derives ultimately 
from the conserved symplectic form in phase space. And candidate states are defined as 
elements of the Fock space built over a one particle Hilbert space whose inner product also 
derives from the conserved symplectic form. 

While we can loosen the form of the equation of motion, it is difficult to know where 
to begin if we do not maintain the existence of a conserved symplectic form. Thus, at least 
for now, we restrict to irregular lattices that nevertheless possess a time layering so that 
the dynamical connection between points proceeds systematically from one time slice to 
the next, allowing a conserved symplectic form to be defined. 

Even when there is a conserved symplectic form, the equation of motion we write will 
generally not admit an initial value formulation, as evolution from initial data may not 
be unique or may not exist in general. This is inevitable when the number of field points 
and/or their connectivity is not conserved from one time to the next, since there can be 
either too few or too many equations to determine the subsequent values. One of the 
interesting questions is just how far such behavior can be tolerated. 

In order to implement the commutation relation and define the Hilbert space in a 
straightforward way we restrict the lattice so that a linear phase space consisting of so- 
lutions to the field equation can be defined. This requires that backwards evolution, if it 
exists, is unique. To include all observables of interest, it is necessary to allow for 'birth 
solutions' that come into existence at a finite time, i.e. that cannot be evolved backwards 
past a certain time. With all of this structure we find it possible to realize the algebraically- 
defined quantum field theory in a way that preserves the equations of motion, has canonical 
commutation relations, and allows for mode birth. The ambiguity of future evolution is 
incorporated naturally into this structure and does not present any fundamental problem. 
We also construct a Fock representation of the algebra, with the help of a particular rule 
for resolving the evolution ambiguity. The Fock representation depends on this rule, but all 
such representations are equivalent, and the theory is fully defined by the algebra without 
the choice of any such rule. The state of the newborn modes is not specified by the earlier 
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data so is extra information not available at an earlier time. 



We study here particle creation in a 1+1 dimensional example, assuming the usual 
vacuum on an early time regular lattice and examining the spectrum of created particles 
in a late time regular lattice after points have been added. Considering all states subject 
to the restriction that there are no incoming particles, we arrive at the conclusion that, 
roughly speaking, one Planck unit of energy is injected per Planck volume (really length 
in the 2d model) created. This calculation serves to illustrate that the theory has been 
defined well enough to carry out such a computation. But it also suggests that we have 
not achieved a viable theory: it would be incompatible with the observed universe to have 
so much energy injected. We return to this issue in the discussion section. 

Other approaches to the problem of formulating mode birth in an expanding universe 



have been explored recently. Kempf |1C] has constructed a theory of fields with finite 
density of states, derived by representing a deformed uncertainty relation between position 
and momentum which implements a minimum spatial uncertainty. A quadratic scalar 
field Lagrangian written with this deformed momentum operator produces a field equation 
whose mode solutions are born at a definite time. The modes cannot be evolved backward 
before the birth time, because the differential equation they satisfy is singular at that time. 
This is closely analogous to what occurs in our discrete model. Unlike our case however, 
forward evolution is not ambiguous. 

In a very different approach, Sasakura [11] formulates discrete field theory of scalars, 
spinors, and gauge fields, on a growing fuzzy 2-sphere. At each discrete time N the scalar 
field is an hermitian operator on the representation space of the spin-iV/2 representation 
of su(2), and the fields 4>n+i, 4>n-, and 4>N-i are related by a discrete analog of the wave 
equation. The classical scalar field equation in an expanding universe with proper time N 
and scale factor \/N + 1 is recovered in a large N, low spin limit. The initial value problem 
is not discussed, but it seems that the field equation does not determine 4>n+i from 4>n and 
<f>N-i- We would guess that there are also birth solutions. The quantum interpretation of 
this theory is not yet entirely clear (at least to us). Both of these approaches share the 
feature that a short distance cutoff exists without sacrificing spatial symmetry, unlike in 
our rather more crude approach. 

The rest of this paper is organized as follows. Section 2 describes the lattice structure 
adopted, and section 3 formulates the action and field equation, as well as the phase 
space description. The ingredients then allow straightforward quantization, as described in 
section 4. In section 5 we work out the example of quantization on a regular 'hyperdiamond' 
lattice, and section 6 deals with the example of a two-dimensional lattice that grows by the 
addition of two points. In particular we construct a Fock space representation and examine 
the resulting particle creation. A discussion of the results, some perspectives, and ideas 
for future work are given in section 7. An appendix contains formulae for the Bogoliubov 

coefficients in the two-dimensional example. We adopt the metric signature (H ) and 

use units with the speed of light c = 1. 
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2. Lattice structure 



We use a type of lattice that is a special kind of causal set (||, EJ. This choice allows the 
discrete structure to carry the metric information in its topology, with no length assign- 
ments. To make the field theory action finite, we restrict to the case where each point p 
has a finite number of causal neighbors, i.e. points q such that there is no point between 
p and q in the sense of the causal order. We denote by IF(jp) the set of future causal 
neighbors, the immediate future. Similarly IP(p) refers to the immediate past. These sets 
are thought of as discrete analogs of a cross section of the future and past light cone in the 
tangent space at a point in the continuum. A pair of causal neighbors is called a link, so 
the structure could be called a finitely linked causal set. 

We choose to preserve symmetry among the future links at a point, in the sense that 
our field theory action is unchanged when they are permuted. In the continuum limit 
this amounts to treating the future links as spacetime vectors of identical length and 
with identical pairwise inner products. Unlike the standard interpretation of the ordering 
relation in causal sets, the links will correspond to spacelike rather than null vectors (except 
in two dimensions), although the convex hull spanned by these vectors will have null faces. 

The action and equation of motion will be well defined on the minimal structure of a 
finitely linked causal set. But to make progress on the quantum field theory, as discussed 
in the introduction, we shall require that the lattice be layered by time slices: the lattice 
will consist of an indexed sequence {S n } of sets of points called slices, such that every point 
lies in exactly one slice S n , and links connect points in slice S n only to points in S n -i or 
S n +l- 

We will use the letters i,j to index elements p 1 G IF(p) and pi £ IP(p) with no 
implication that the number of elements in any of these sets is equal. Multiple indices are 
also allowed, thus for example 

pj := ( Pi ) j G IF(IP(p)). (2.1) 

Note that the index order matters: pP i G IP(IF(p)) is generally not equal to p^ . One 
calls p l a child and pi a parent of p. Then p^ is a sibling while p?i is a mate. All manner 
of polygamy and incest is permitted by the layered structure, provided it is not cross- 
generational. 

To carry out the quantization in close analogy with linear quantum fields in the con- 
tinuum, we will make one further restriction on the lattice structure: backward-in-time 
evolution using the field equation we adopt, if it exists, is unique. 

2.1 Hyperdiamond lattice 

The regular prototype of the lattices we will consider is a d-dimensional hypercubical space- 
time lattice tilted on its vertex so that a preferred time-like direction runs from one corner 
of the cube to the complete opposite. We call this the hyperdiamond lattice after the 

1 In the causal set literature a "slice" refers to a maximal set of unrelated points. Our slices here consist 
of unrelated points, but they need not be maximal. For example the points e and / in Fig. |l| are unrelated 
to points on the previous slices, so the latter are not maximal. 
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terminology of Ref. [12, p^| . A step along one lattice link runs along one of the cube's 
edges, and represents a movement in the preferred time-like direction and a movement in 
the orthogonal space to one of the d vertices of a [d — l)-dimensional tetrahedron centered 
on the original point. The lattice thus possesses a tetrahedral rotation symmetry. 

We choose the step-speed so that the future continuum light cone at a point lies inside 
and tangent to the polyhedron formed by the set of future links. With this choice the 
hyperplane faces of the polyhedral cone are null hyperplanes and the links are spacelike. 
This determines the step speed as follows. 

The d link vectors are given by 

n a i= t a + as* (2.2) 

where t a is a unit timelike vector along the hypercube diagonal and, for each i, sf is a 
spatial unit vector orthogonal to t a pointing to one of the vertices of a (d — l)-dimensional 
tetrahedron. By tetrahedral symmetry, the sum ^ i s" must vanish. If the cube faces are 
to be null, then their normal vector must lie within them and must be a null vector. Again 
by tetrahedral symmetry this normal for the face omitting nf must be the sum YLj^i 
of the d — 1 link vectors defining the face. Since the d spatial vectors sum to zero this is 
equal to (d — l)t a — asf, which is null provided a = d — 1. Hence the link step speed is 
d — 1, which is the speed of light only in two spacetime dimensions. 

Squaring the relation Yli s i = reveals the spacetime inner products of the spatial 
vectors, sfsj a = — 1) for i ^ j. With our signature, the spatial inner products have 
the opposite sign, so the cosine of the angle between any two is —l/(d— 1). The inverse 
spacetime metric can be expressed in terms of the link vectors as 

An equivalent description of the hyperdiamond lattice is obtained beginning with a 
system of d null coordinates v l , i.e. coordinates whose level sets are null hyperplanes, such 
that the metric is invariant under permutations of these coordinates. The corresponding 
basis one- forms dv l are null, and the dual basis vectors (d/dv l ) a are equal to the nf defined 



in (2.2), when the v l are scaled so that t = ^ v l is a proper time coordinate. Discretizing 
the v % then yields the lattice. 

The d-dimensional hyperdiamond lattice can be compactified to have spatial sections 
that are (d — l)-tori. Consider for illustration the case d = 4. Then the immediate future 
of any point is a tetrahedron. Translation along any edge of the tetrahedron is clearly 
a symmetry of the hyperdiamond lattice, since every point is equivalent. We can thus 
identify the points related by any integer number of such translations along three edges at 
one vertex to obtain a 3-torus. One can also compactify to "twisted tori", or even Klein 
bottles. 

3. Lattice dynamics 

A continuum massless, minimally coupled scalar field in d spacetime dimensions has the 
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action 



S = \ J d d x^-gg ah d a ^d h 4>. (3.1) 



We first discretize this on a hyperdiamond lattice, then use the result to motivate an action 
on an irregular lattice. After discussing the resulting equations of motion, we introduce a 
phase space and symplectic form to be used in the quantization. 

3.0.1 Hyperdiamond lattice action 

Given a hyperdiamond lattice with edge links enf scaled by a small length e, we insert the 



metric decomposition ( |2.3| ) in the action (3.1) and replace the partial derivatives by finite 
differences 

en?0 o 0->0(p + eni)-0(p). (3.2) 

This yields the discrete action 

5 reg = V^e d ~ 2 \ did^T) {^ pi) ~ ^ p) ) " m ) ' (3 ' 3) 

where V^e d is the spacetime volume of one hypercubical cell of the lattice, and p l denotes 
the point future- linked to p in the direction n™. 

One could instead use the past links rather than the future links, or even a combination 
of the two, in the discrete action. We use the future links for no good reason. 

3.0.2 Irregular lattice action 

Using the regular case (|3.3|) as a guide, we posit an action for a scalar field on a general 
causal lattice. We cannot copy that exactly at each point p, since the number of points 
dp in the immediate future IF{p) varies from point to point, hence different powers of 
the dimensionful quantity e appear. Our definition of the discrete theory simply ignores 
the powers of e and the volume factor that goes along with them. We preserve the 
denominator d p (d p — 1), which corresponds to the number of terms in the sum, so that the 
action at p is always the average of the products of the differences along the future links 
at p. With these choices, the irregular lattice action we adopt is 

= 2 E d u- 1) ^ ~ m ) ~ m) ■ (3-4) 

i¥=3 

3.1 Discrete field equations 



Varying the field in the irregular lattice action (3.4) one obtains the discrete field equation 

U(j) = 0, (3.5) 
where the discrete D'Alembertian operator □ is defined by 

□<Kp) = -Hp) + E T ^ + E Hpi) ~ E a .(d.-i) HpiJ) - (3 ' 6) 

i p i Pi ij Pi\ Pi / 



-7- 



In the regular case, the D'Alembertian is the average over the children, plus the average 
over the parents, minus the average over the siblings, minus the value of the field at 
p. It does not have this interpretation in the general case, because the weights of the 
parent and sibling terms are different. It can be rewritten to look more like the continuum 
D'Alembertian df — V 2 : 

E I - *oo) - E £ {m - «*)) - E s^tt) M - *)) • <"> 

i 1 i yi i,j rt\ ft i 

In the regular case the first two sums are the average future and past time-differences, so 
their difference gives a discrete version of the second time derivative. The third sum is the 
average of the spatial neighbors minus the value at p, which is a discrete version of the 
spatial Laplacian. 



3.2 Initial value problem 

The field equation (3.5,3.6) involves values of the field on three consecutive time slices. The 
equation that arises from varying (j){p) on the middle slice determines the average value on 
IF(p) given the values on the first two slices. If the number of points on the middle and 
third slice is A2 and A3 respectively, then one has A2 linear equations in the A3 unknown 
values on the third slice. Supposing N2 and .A3 are finite, the future solution is not uniquely 
determined if A3 > A2, and there is generically no future solution if A3 < A2. Existence 
of a future solution when A3 > A2 depends upon the connectivity of the lattice. 

An example we shall make use of is a two-dimensional regular diamond lattice, where 
each point has two parents and two children, and periodic boundary conditions are imposed 
so that at each time there is a total of A points arranged on a circle. The field equation 
determines the sum of the values on each adjacent pair of points, so if A is even (or infinite) 
the evolution is determined only up to the addition of a spatially-alternating function with 
arbitrary time-dependence. If A is odd the evolution is uniquely determined since the peri- 
odicity is inconsistent with the alternation. It is possible to regard the alternating functions 
in the even or infinite cases as "pure gauge", and to formulate a corresponding dynamics 
of the gauge-invariant degrees of freedom, but we shall avoid this extra complication by 
requiring A to be odd in the two-dimensional example we study. 

Note that even if the number of points on each slice is the same, the equations could be 
over-determined so that no future solution would exist. For example this could happen if 
the future of some set of points on one slice contains fewer points on the next slice. But even 
if these numbers match there could still be trouble. Consider the two-dimensional example 
just discussed in the case that A is even. The equations take the form <f>(i) + <j>{i + 1) = Hi, 
with i = 1, . . . , A and 4>(N + 1) = (f)(1). The sum of the left hand sides of the odd-z 
equations is equal to the sum of the even-i ones. (For example, with A = 4 one has 
[<p(l) +<p(2)] + [(f)(3) +0(4)] = [(f)(2) + (f)(3)] + [(f)(4) + (f)(1)].) This implies that there can be 
a solution only if the condition Y^ odd = Seven H% holds. It happens that this condition 
is automatic if the lattice is a regular diamond lattice, but it could fail to hold on a lattice 
with altered connectivity prior to the slice in question. 
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3.3 Phase space of half-solutions 

Since initial data cannot in general be evolved indefinitely or uniquely, it is not convenient 
to define a phase space in terms of initial data. Instead one can try to use the covariant 
approach that identifies phase space as the space of solutions. 2 This takes care of non- 
uniqueness of evolution, as two distinct evolutions simply correspond to different points in 
phase space. Non-existence of evolution is another matter. 

One might try to associate to each region a different phase space, consisting of the 
solutions that are defined everywhere in that region. This idea may be compatible with 
that of a net of local algebras on which the algebraic formulation of QFT is based ||. 
Another approach might be to introduce just one phase space, without a linear structure. 
But, for this initial investigation, we choose to exploit the linear structure of a single phase 
space in strict analogy with standard linear quantum field theory. Thus we seek a definition 
of phase space as a single vector space. In particular, we require an operation of addition 
of two solutions that defines an abelian group, i.e. that is commutative, associative, and 
possesses an identity and an inverse. This requirement leads to further demands on the 
partial solutions and lattice structure. 

Let </>(s..t) denote a solution that exists between the discrete times s and t and is 
undefined outside that interval. The additive inverse of cj){s..t) must surely be —cj){s..t), 
but these two functions are defined only on the interval [s, t], whereas the identity must be 
the zero solution which exists for all time. We thus learn that the operation of addition 
must include a maximal extrapolation (using the field equation) beyond the original range 
of definition of the sum. However, since evolution is generally not unique, extrapolation is 
not well-defined. To cope with this we restrict to lattices for which backward evolution, if it 
exists, is unique, and we define the phase space to consist of solutions that are defined for 
all time, together with solutions that are defined on a semi-infinite interval 3 [i, oo). These 
latter we call half- solutions with birth time t. On this space, the operation of "sum on the 
common domain, then maximally extend to the past" defines an abelian group. To see 
that this is associative consider (a + b) + c and a + (b + c). In the region where a, b, and 
c are all individually defined these sums are both equal to the standard sum a + b + c. 
Globally, they are both the maximal past extension of this sum, hence they are equal. 

As noted in section ( |3.2[) , there is no guarantee that all — or even any — initial data can 
be extended indefinitely into the future. Despite this we can define the phase space using 
the half-solutions. If there are none, then the theory is empty. 

3.4 Symplectic form 

In a general Hamiltonian context, a symplectic form f2 on phase space is a non-degenerate 
closed two-form ||. For a linear phase space T the tangent space at a point is naturally 
identified with the phase space itself, so can be treated as a map from T x T to the real 
numbers. In the present application, Q maps pairs of half-solutions to numbers. 

2 See H for an introduction to the covariant phase space formalism and its application to quantization 
of linear field theory. 

3 The theory could also be set up using solutions defined on intervals of the form [t, to], where to is a 
fixed finite time. 



- 9 - 



To motivate the discrete version, recall first how it works in the continuum. For a 
real Klein-Gordon field in D-dimensional space-time, the symplectic product between two 
solutions ijj and 4> is defined by 

n (v, 4)= [u d a <t> - <rv) ds a , (3.8) 

Js 

where S is a Cauchy surface. This product is well-defined on the covariant phase space since 
it takes the same value when evaluated on any surface S provided the field equations are 
satisfied. A simple way to verify this fact is to integrate the expression ipDcp — <f>nif), which 
vanishes by virtue of the field equation, over a space-time region between two Cauchy 
surfaces. This integrand is the total divergence ^ a {ipd a 4> — 4>d a iJ)), hence the volume 
integral is equal to the difference of ( |3.8P between the final and initial surfaces (assuming 
no contribution from a spatial boundary). Since the volume integral vanishes, the two 
spatial integrals must be equal. 

The same construction works in the discrete case, with □ replaced by the discrete 
version ( |3.6| ). The sum of (pHip — ifiHcj) over a sequence of slices reduces by cancellations to 
a pair of boundary contributions on the initial and final slices (assuming no contribution 
from a spatial boundary). This shows that the discrete symplectic form 

nfcM) = x(^frW) - HpW)) (3.9) 

takes the same value when computed on any slice S in the domain of definition of two 
solutions tp and <p to the discrete field equation. It is thus well-defined on phase space. On 
a spatially infinite lattice some kind of falloff condition on the fields would be required to 
ensure that the infinite sum defining the symplectic form converges. 

3.4.1 Degeneracy 

The symplectic form for the continuum Klein-Gordon field is non-degenerate, i.e. £l(ip, •) = 
only if ip = 0. This is not necessarily true on the lattice. For example, suppose A is any 
spatially alternating configuration on a two-dimensional periodic even diamond lattice. 
Then, as discussed in section |3.2| , A satisfies the field equation. It is thus an element 
of the phase space, so Q(A, •) is well defined, and is in fact equal to zero. Whenever a 
symplectic form has a vector field Y for which £l(Y, •) = 0, the Hamiltonian vector field 
Xh is determined by Hamilton's equations dH = Q(-,Xh) only up to addition of an 
arbitrary multiple of Y. In the present case, this corresponds to the statement that to 
any solution eft one can add an arbitrary spatially-alternating configuration to obtain a 
new solution with the same initial data. This indeterminacy is a gauge freedom, and the 
physical phase space is the quotient of the solution space modulo the gauge freedom. 

3.5 Poisson brackets 

The symplectic form, if not degenerate, provides the definition of Poisson brackets between 
any two phase space observables, i.e. functions on phase space. Given two such functions 
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/ and g, the Poisson bracket is a third function defined by 



{f,g} = n ab d a fd b g. (3.10) 

Here we use (abstract) index notation, and Q, ab is the inverse of the symplectic form flab, 
defined by fl ab fl bc = 5 a c . 

In the linear case, to each phase space point <j> is associated the linear observable whose 
value at tp is fl((f),ip) = fl ma (p m ip a . This function has the gradient d a fl(<p, •) = fi ma m ; 
therefore 

{n(^o, n(0,-)} = -n(^). (3.11) 

It is this relation that we quantize. If fl is degenerate, an algebra of gauge-invariant 
observables can be defined directly by the relation ( 3.1 1| ) , which does not involve the 
inverse of fl. 

3.5.1 Local fields and locality 

Though not central to our study, we shall make a few remarks here about local observables 
and their Poisson brackets. This section can be omitted without missing anything essential 
to the central focus of our paper. 

First, note that evaluation of solutions at a lattice point x defines a linear map <p — > 
4>{x) only on those half-solutions whose birth time lies before x. Since this map is not 
defined on the entire phase space it is not strictly speaking an observable. The Poisson 
bracket {(f)(x), (j>(y)} is thus not globally defined on the phase space, but it should be 
meaningful locally in phase space. 

To discuss the locality in space-time of the Poisson brackets, let us restrict to the case 
of a compact (i-dimensional hyperdiamond lattice, chosen so that evolution from initial 
data always exists and is unique. Then the local field 4>(x) is an observable. We call this 
a deterministic lattice. 

In the continuum, the Poisson bracket of two local fields is given by 

{<j>{x), 0(y)} = G R (x, y) - G A (x, y), (3.12) 

where Gr and Ga are the retarded and advanced Green's functions of the Klein-Gordon 
equation. If x and y are spacelike related this bracket vanishes. On the lattice (in greater 
than two dimensions) the links are spacelike, so we do not expect the bracket to vanish at 
all spacelike separations. However, not even a lattice locality property survives, as we now 
explain. 

On a deterministic lattice fl is necessarily invertible, 4 so there exists a solution F x such 

that 

<f>(x) = n(F x ,ct>). (3.13) 



4 If A is a solution such that tl(A, ■) = 0, then the future and past averages of A must vanish at every 



point. The field equation (3.E,3.6) then implies that A = 0, unless d = 2. In the case d = 2, uniqueness 
of evolution requires an odd number of points on each slice, which together with the vanishing of averages 
implies A = 0. 
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Evaluating on the slice S x , one sees that in order for Eq. ( 3.13 ) to hold for all 4>, the 
solution F x must vanish everywhere on the slice S x containing x, have a mean value of — 1 
on IF(x), and have a mean value of zero on IF(v) for every other point v in the slice S x . 
The bracket of local fields is thus 

y>{x),^{y)} = -n(F x ,F y ) = -F y {x). (3.14) 

The second equality follows immediately from the defining property ( |3.13j ) of F x . If x 
and y lie on the same slice then F y {x) = 0, so the equal-time bracket vanishes, as in the 
continuum. If y lies on the subsequent slice, we need to use the field equation to evolve 
F y back one step to the slice S x in order to evaluate F y (x). The field equation ( |3.5| , j3~li| ) 
on S y implies that the average of F y on IP(y) is +1 and the average on IP{w) vanishes 
for every other point w in the slice S y . On a hyperdiamond lattice this clearly forces F y 
to have support on S x arbitrarily far from y, hence there will be many such pairs x, y for 
which Fy(x) ^ 0. 



This nonlocality of the bracket ( 3.14 ) survives quantization, so it implies arbitrarily 
spacelike distant observables fail to commute. It is not clear that this is really a problem, 
however, since the average of field values on each IF{x) may have purely local brackets. 
This average corresponds to the observable £l(6 x , •), where 5 X is equal to 1 at x and vanishes 
everywhere else on the slice S x containing x and the subsequent slice. That is, we have 



(x l )=n(5 x ,<P), (3.15) 



so the bracket of two such local averages is 

{ \ £ \ £ W)} = S v) = -\Y, 5 v^)- ( 3 - 16 ) 



Unlike for F x , the data defining 5 X is localized. Since the field equation (3^,3.6) is local 



this means that the support of the solution S x might propagate "causally" away from x 
to the future and the past. If so, then the bracket of the averages is non-zero only if x 
and y are causally related in the sense determined by the local field equation. On an odd 
two-dimensional lattice the propagation of 5 X is indeed causal. We do not know if that 
property holds in higher dimensions. 

4. Quantization 

To quantize one replaces the classical observable by an element of a C*- 

algebra possessing the canonical commutation relations (CCR), 

n ((/), •) , n (v>, ■)] = (<t>, V>) i, (4- 1 ) 

and such that ■)* = (Hereafter we set h = 1.) The quantum theory is then 

completed with the choice of a state, i.e. a linear functional ui on this algebra satisfying the 
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positivity condition co(a*a) > 0, with a an arbitrary algebra element, and the normalization 
condition uj(1) = 1. This notion of state adapts the concept of a density matrix to the 
most general setting. 

Here we follow the simple conventional procedure of constructing the class of states by 
representing the algebra as operators on a Fock space. The Fock space is built over a "one- 
particle" Hilbert space which is defined using the Klein-Gordon inner product between 
complex solutions to the classical field equation: 

<v, 0) = <f>) = «v ±- {4>{pW) - mW)) ■ (4-2) 

P es-i ap 



The definition ( |3.9| ) of Q is extended here by complex linearity to the space of complex 
solutions to the field equation. Although this is not a positive-definite inner product it can 
be used to construct the Hilbert space as follows. 

Choose a positive-norm subspace which is orthogonal to its conjugate, and which 
together with its conjugate spans the space of solutions. Let be an orthonormal 

basis adapted to this decomposition, so that 

= Sij = - (&,&>, (d, & = 0, (4.3) 
and define the field operator ij> as 

4> = J2(& a i + & a i): (4-4) 

i 

where 

[ai,aj] = 5ij, [ai,aj] = = [at, a]]. (4.5) 



Then the commutation relations ( |4.l|) are satisfied if the quantized linear observables are 
represented by 

[n(V, •)] rep = inw, 0) = » J] ( n M &) a * + n M &)<4) • ( 4 - 6 ) 

i 

The 'Fock vacuum' state |0) is defined by the condition Oj|0) = for all i, and the Fock 
space is the Hilbert space spanned by all finite norm vectors obtained by acting on |0) with 
products of the raising operators at. The field operator (4.4) is self-adjoint acting on the 
Fock space, as befits a real field. The state vectors of the theory are the elements of Fock 
space, while the algebraic states are density matrices constructed therefrom. 

The construction in the last paragraph depends upon the basis but on a spa- 

tially finite lattice the representation of the algebra is independent of this choice. This 
follows from the Stone- Von-Neumann theorem, which states that all representations of 
the canonical commutation relations for a finite dimensional phase space are unitarily 
equivalent |9|. 
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5. QFT on a hyperdiamond lattice 



In this section we work out the example of a toroidally compactified hyperdiamond lattice, 
beginning with the classical fields and then quantizing. We initially consider an arbitrary 
dimension d. We find that the discrete field equation has exponentially growing modes 
except in two dimensions. In the next section we use these results in the two-dimensional 
case to determine the particle creation on a growing lattice. 

5.1 Hyperdiamond modes 

We can use the harmonic modes 

Vk=N k e- lk * v \ (5.1) 



where v % are the null coordinates discussed at the end of section 2.1, and Mk is a normaliza- 
tion factor. We adopt units with the lattice spacing e = 1, so one future link on the lattice 
corresponds to adding unity to one of the d null coordinates v % . The toroidal boundary 
conditions select a set of allowed wave vectors h. 



Inserting into the field equation (|3.5|,p7q) yields the dispersion relation 



' ' x ' 7 ' did 1) 

For small ki (compared to 1/e) this reduces to ^2i^j hkj = 0, which says that ki is a null 
vector. That is, the continuum dispersion relation is recovered in this limit. The special 
cases where all but one of the ki vanish are exact solutions both on the lattice and in the 
continuum. These are plane waves whose constant phase surfaces are level sets of one of 
the d null coordinates. 

If the mode is written as exp(— ioot + ik ■ x), the frequency u and spatial wave vector 
k are related to ki via 

ki = uj — (d — l)k ■ §i, (5-3) 

where the unit vectors Si point to the corners of the spatial tetrahedron (see discussion in 
section |2.1| ). We do not assume u is real since we wish to discover whether or not there 
are modes that grow or decay exponentially with time. However, as in the continuum, we 
restrict to real spatial wavevectors. 

For general ki the dispersion relation ( |5.2| ) can be expressed as 

SW+ E(-*, = «^I±^1 (5.4) 



where 



In terms of frequency and spatial wavevector X can be factored as 

E(u,k) = e iuJ z(k), (5.6) 
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where 

z (k) = I^ e -i(rf-l)^\ (5.7) 

i 

Using the polar form z(k) = re tS , the dispersion relation (|5.4|) becomes 



dr 2 + d - 2 . , 

Since the spatial vector is assumed real, < r < 1. 

In the two-dimensional case (5^) reduces to cos a; = cos A;, where k now refers to the 
spatial wavevector. If cos k ^ the general solution is uj = ±k, as in the continuum, 
corresponding to right and left moving modes. Equivalently, &2 = or k\ = 0. If cos k = 0, 
then the division by r in going from ( |5.4|) to ( |5.8| ) is disallowed, and in fact k = ±7r/2 is also 
a solution for any value of u. These are spatially alternating modes, since the spatial step 
length is 2e. They are the pure gauge solutions discussed in Sec. 3^, which are excluded 
on a periodic lattice with an odd number of points. 

If the dimension d is greater than two, then the right hand side of ( |5.8| ) is a positive 
function that is concave upward and diverges as r goes to zero or infinity. This function 
is less than unity only for (d — 2)/d < r < 1. When the wavevector is small, compared to 
the inverse lattice spacing, r is close to unity, so the frequency is real. If it is possible to 
choose k so that r < (d — 2)/d, then there will exist a complex solution to (|5.S|) , implying 
the existence of exponentially growing or decaying modes with wavelengths of order the 
lattice spacing. We have not tried to determine whether there exist boundary conditions 
that exclude such modes, but we now show by an example that they can exist. 

To simplify the computations we choose the spatial wavevector to be orthogonal to 
§i for i > 3, so it has the form a(§i — §2) for some constant a. We then find z(k) = 
[(d — 2)/d] + (2/d) cos ad. Hence r = \z\ < (d — 2)/d if cos ad is negative. For such values 
of a, the mode is unstable. 

On a compact lattice, such a mode must satisfy periodic boundary conditions. Consider 
a toroidal compactification of the lattice as described at the end of Sec. ^l], with identifica- 
tion after N lattice translations in each of the (d — 1) directions Sj — §1 for i = 2, . . . d. The 
lattice step with i = 2 for example is given by (v 1 , v 2 ,v 3 , . . . , v d ) — ► (v 1 — 1, v 2 +l, v 3 , . . . , v d ). 
In order for the mode function to be periodic after N steps in each of these directions, its 
wavevector must satisfy k% — k\ = (d— l)k ■ (si — Sj) = 2irrii/N. With the above form of k, 
these periodicity conditions reduce to 2ad = 2irri2/N, and ad = 2irni/N for i > 2. Thus 
we must choose the rtj equal to a common value n for i > 2, and 122 = 2n. This mode 
then satisfies the periodicity condition if ad = 2irn/N. If we choose n so that cos2im/N 
is negative, this mode is unstable. 

We are not interested in trying to make sense of the quantum theory in the presence 
of unstable modes. Perhaps under appropriate boundary conditions they are excluded. 
Certainly they can be eliminated by adopting a different sort of discretization scheme. 
The problem need not concern us here, since the example we shall study in detail is two- 
dimensional. Moreover, our general approach to quantizing the field in the presence of 
mode birth is independent of the issue of stability of the chosen finite difference scheme. 
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In the continuum, for a fixed spatial wavevector there are two frequencies that solve 
the dispersion relation, to = ±\k\. On the hyperdiamond lattice, there are still two fre- 
quencies (modulo 2ir) that solve the dispersion relation, but, except in two dimensions, the 
frequencies are not negatives of each other. Indeed, if (oo,k) is a solution to ( |5.8| ), then so 
is (a/, k) if 

J = -to -26. (5.9) 

This is equivalent to the condition £(u/, k) = £(— u, —k) = £(w, k), where the last equality 
holds only if to is real. 

Another difference with the continuum in more than two dimensions is that if (to, k) is 
a solution, then (uo, —k) is generally not a solution. Under k — ► —k we have z(k) — ► z(k), 
so (r, 6) — > (r, —0). Therefore — A;) is a solution if and only if (uo — 8) = ±(u> + 6), i.e. if 
6> = or if to = 0. 

As in the continuum, a massless scalar field has a pair of spatially constant solutions 
on the lattice. To complete the list of modes we must include these zero-modes. The 
constant field with ki = is a zero-mode, as is the mode that grows linearly in time, 
<f) oc v % . These play a role in the particle creation calculation carried out below for the 
two-dimensional growing lattice. 

5.2 Quantization on a hyperdiamond lattice 

To carry out the Fock space quantization procedure we must identify the decomposition of 
the space of complex classical solutions into positive and negative norm subspaces. This 
yields the representation of the algebra, and a pure state is specified using an element of 
the Fock space. 

5.2.1 Harmonic modes 

The Klein-Gordon norm (4.2) of the modes (5.1) is 



(rj k ,r, k ) =2N\jV k \ 2 lmE, (5.10) 

where ./V is the number of spatial points on the lattice. We assume henceforth that all 
components of the wavevector ki are real, so that fj k rj k = |A4| 2 - Modes with ImS > have 
positive norm, and those with ImS < have negative norm. The normalization factor (up 
to an arbitrary phase) is 

Afk =(2iV|rm£|)~ 1/2 . (5.11) 

The space of modes thus decomposes into these two subspaces, which are conjugate to each 
other. Moreover, they are mutually orthogonal. Consider the pairwise KG inner products: 

( Vk , Vl ) = i(E(i) - E(fc)) J^VkVi- (5-12) 

pes 

The sum will vanish unless the spatial wave vectors for ki and li are equal. As seen in 



Sec. 5.1, for a given spatial wavevector the allowed frequencies are uo and — oo — 29, corre- 



sponding to S(/c) and S(/c') = £(&;). Thus if the modes with the same spatial wavevector 



- 16 - 



are different, then = S(Z). In this case, one mode has positive norm and the other 
has negative norm, and the inner product vanishes. The modes thus provide a suitable 
orthonormal basis of complex solutions for constructing a Fock space representation of the 
field algebra. We use the notation for the positive- norm modes, which are just those % 
for which ImS(fcj) > 0. Explicitly, they are given by 

V2iVImS(A; i ) ' K ' ' 

5.2.2 Zero-modes 

The constant solution p = 1 and the linear-in-time solution a = t = ^ ■ v l both have 
vanishing norm, and they are orthogonal to all the harmonic modes % for k ^ 0. The 
KG inner product between them is (p, a) = iN . This two-dimensional zero-mode subspace 
can be split into a positive-norm direction and its conjugate, orthogonal, negative-norm 
direction, but there is no unique way to do so. The different choices of positive-norm 
zero-modes are p — ijcr, where 7 = 7i? + iyi is a complex number with positive real part. 
Explicitly, the normalized, positive-norm zero-mode choices are 

5.2.3 Field operator and Fock vacuum 

The field operator is given as in (|4.4j) by 



4> = ^2{£kak + £fc4)> ( 5 - 15 ) 

k 

where the sum is over those ki that satisfy the dispersion relation with Im£(/cj) > 0, and 
the 0^ and at satisfy the usual commutation relations of lowering and raising operators. 

The Fock vacuum is the state satisfying ak\0) for all k. In this state, the expectation 
value of the field operator vanishes, and that of the square yields 

(O|0(p)0(p)|O) = (2iV 7K )- 1 ((l+ 7/ t) 2 + ( 7i? t) 2 ) +^r(2ATm£(£:)) _1 . (5.16) 

The time dependence coming from the zero-mode shows that the state is not time inde- 
pendent except for 7 = 0, which does not correspond to a normalizable zero-mode. The 
field amplitude of the zero- mode is like a free particle, whose wavefunction always spreads. 

One might be inclined to identify the state with the lowest energy as the preferred 
vacuum. As for a free particle, there is no normalizable state that corresponds to the 
lowest energy of the zero-mode degree of freedom. Considering just the other modes, there 
is also a fundamental issue on the lattice: since time evolution is discrete, there is no 
uniquely defined Hamiltonian, but rather a finite-time evolution operator U(At). Taking 
the logarithm one could define a Hamiltonian H = (i/At) InC/, but this is ambiguous due 
to the freedom to add a multiple of 2ir/At to any eigenvalue without changing the evolution 
operator. 
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Even if there is no meaningful lowest energy state, one can still ask if the Fock vacuum 
is an eigenstate of the evolution operator. Other than for the zero-mode, the answer is 
yes: if the evolution operator U is defined by W^v 1 , v 2 , . . . , v d )U = 4>{v l + 1, v 2 , . . . , v d ) 
then, omitting the zero-mode, U\0) = |0) (up to an arbitrary phase). Combining such 
translations of the null coordinates one obtains all space and time translations. 

6. Mode birth on a growing two-dimensional lattice 

To explore the phenomenon of mode birth in a simple setting we now specialize to two- 
dimensional lattices. We consider a model in which the lattice is initially regular and then 
some points are added, after which the lattice is again regular. We first specialize the 
results of the previous section to two-dimensional regular lattices, then turn to the case 
where points are added. 

6.1 Regular diamond lattice 



On a regular diamond lattice in two dimensions the field equation (3.5,376|) can be written 
as 

(<p( P L ) + <Xpl)) - {Hpl l ) + Hp))} + [{Hp r ) + Hpr)) - (<f>(p) + Hpr r ))} = o, (6.1) 

where p and p R are the left and right children, pl and pr are the left and right parents, 
and p L L and p R R are the left and right siblings of p. In terms of the diamond operator, 

<>Hp) = (4>(p l ) + Hpl)) - {Hp l l ) + Hp))), (6.2) 

the field equation (|6.1| ) takes the form 

OHP) + <>HPR R ) = 0. (6.3) 

This asserts that the sum of every pair of adjacent diamonds vanishes, so the diamonds on 
each time slice have the same magnitude and alternating sign. On a finite periodic lattice 
with N points there are N diamonds. If N is odd the alternation is consistent only if each 
diamond vanishes individually. 



The mode functions were found in Sec. 54. For the right-moving modes (&2 = 0) we 



have 2ImS = smk\. The positive-norm modes ( 5.15 ) are those for which sinfci > 0, and 
are given explicitly by 

6c = J^ e-^ 1 . (6.4) 
V-/V sin k± 

A similar expression with (k±, v 1 ) replaced by (k2,v 2 ) yields the left-movers. The zero-mode 
(5.14) takes the same form in all dimensions. 



As explained in section 3.2, we restrict to odd values of iV to avoid the alternating 
gauge freedom. Periodicity for the mode solutions ( |6.4| ) implies that k% = 2irn/N for some 
integer n, so the spectrum of modes is discrete. We obtain the complete set if n runs from 
— (N — 1)/2 to (N — 1)/2. The positive-norm modes are those with positive k\ in this range, 
and the time-independent zero-mode corresponds to n = 0. There are N positive-norm 
modes all together: [N — l)/2 left-movers, (N — l)/2 right-movers, and a zero-mode. 
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6.2 Growing lattice: mode birth 

We now consider a lattice which is irregular in a localized region where some points are 
added. To preserve an odd number of points an even number must be added on each slice. 
We add just two points, with the connectivity 5 illustrated in Fig. [l|. 

The point labels a, b, c, d, e, f, g, h will be re- 
ferred to in this section. The equation of motion 
from varying the field at the birth point e implies 
that 

0(e) = (0(a) + 0(6)) /2, (6.5) 

i.e. 0(e) is equal to the average of the future linked 
values. A similar equation holds for the birth point 
Figure 1: Structure of the lattice /• 

around the birth event. Evolution of solutions backward through the 

birth event, if possible, is unique, as required by 
our linear phase space construction described in section ^D]. Evolution forward is always 
possible, but is not unique since the values at the newborn points are undetermined. All 
such evolutions correspond to points in the phase space. 

We call the slice containing the birth points the birth time. Evolving backwards, the 
values at the birth points can be determined from those on the subsequent two slices, but 
unless it so happens that the birth point values are the average of their children's values, 
the field equation at the birth points is not satisfied. Put differently, there may be no way 
to take the next step backwards that satisfies the field equation at the birth time. Hence, 
there are half-solutions (cf. section |3.3| ) . The initial data for these birth solutions are given 
at the birth time and on the subsequent slice. 

The general definition of the quantum theory in section || applies to this growing 
lattice. The phase space includes the whole solutions together with the half-solutions. 
The algebra of linear observables on this phase space is canonically quantized, yielding the 
quantum field algebra. Together with a choice of state — i.e. a positive linear functional on 
this algebra — the definition of the quantum theory is complete. To construct a concrete 
representation, define a particular state, and compute the value of some observables, we 
now use the method of mode decompositions. 

6.3 Mode decompositions and Fock space 

After the birth event the lattice is regular, and therefore the Fock space can be based on 
the harmonic modes (|6-4 ) together with the zero- mode ( |5.14| ). A simple set of observables 
to examine are the number operators Nk = atafe for these 'out' modes. We shall compute 
the expectation values (0|iVfc|0), taking for the state |0) the Fock vacuum based on a set of 
'in' modes. 

To find a complete set of in-modes we begin with the modes before the birth event, and 

choose a particular evolution among the many that carry these modes through the birth 

event. We then determine the half-s olutions that are orthogonal to the evolved in-modes, 
5 We also investigated the case where e and / are children of h, which yielded similar results. It might 
be interesting to explore the consequences of more complicated connectivity of the birth points. 
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and augment the Fock space with these. This requires only the values of all modes at the 
birth time and the subsequent slice. 

It is worth emphasizing that the arbitrary choice resolving the ambiguity of evolution 
through the birth event does not affect the quantum algebra of observables, since that 
is constructed from all solutions and half-solutions. Nevertheless, the evolution rule does 
affect the inner product of the evolved pre-birth-modes with the half-solutions, since the 
inner products must be evaluated on or after the birth slice. One cannot evaluate on a 
slice prior to the application of the arbitrary rule, since the birth solutions are undefined 
there. 

Although it does not affect the algebra as a whole, the evolution rule does have a 
physical impact in our construction through the definition of the Fock vacuum. The Fock 
space is built from an orthonormal basis of positive-norm modes that are orthogonal to 
their conjugates. The evolution rule affects the construction of this basis, so it affects the 
Fock vacuum. Since all representations are unitarily equivalent however, any Fock state 
defined using one evolution rule can be expressed as a Fock state in a different Fock space 
constructed with a different evolution rule. Thus we are not choosing between theories, 
but rather just between equivalent representations. 

6.3.1 Evolution rule 

For simplicity here we choose the evolution rule that sets the evolved field to zero at the 
birth points e and /. It might seem more natural to linearly interpolate the field between 
the neighbors d and g of the birth points. However, as just explained, this would change 
only the Fock vacuum, not the theory itself. We did examine the particle production using 
the Fock vacuum arising from the "interpolation rule" , and the results were similar to those 
of the "zero rule" . Since the details are simpler for the zero rule, we adopt it for expository 
purposes. 

Let us now see how the fields evolve through the birth event using this evolution rule. 
With cj)(e) = 0, the field equation at the birth point e implies that <ft(a) = —<p(b). 

Similarly the equation at / implies that 0(c) = —<p(b). Taken together these imply that 
0(a) = 0(c). 

Now consider the rest of the equations at the birth time. The field equation at the 
neighbors d and g differs from the standard form (|6.3| ) only in the shape of one of the two 
diamonds. In particular, in the equation from the variation of the value 0(d), the diamond 
on the right consists of the vertices (d,a,g,h), while for the variation of <p(g) the diamond 
on the left consists of (d,c,g,h). These two diamonds agree except at the top vertices a 
and c. The diamonds alternate sign around the compact lattice, and there are an even 
number N + I of diamonds (since in effect the middle diamond has split into two, hence 
the two middle diamonds must be opposite). But we already saw that the field equations 
at the birth points imply 0(a) = 0(c), so these diamonds must agree. The only consistent 
solution is that all the diamonds vanish. With this information we can now easily evolve 
the pre-birth harmonic modes and zero-modes to the next slice after the birth time. 
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6.3.2 Choice of birth modes 

There are two birth points, so there are four independent real birth modes. We wish to 
choose these so they are orthogonal to the evolved pre-birth modes. The birth modes are 
to be characterized by their initial data at the birth time, i.e. their values at the birth time 
and on the next slice. 

To identify these birth modes, we first note that any mode whose initial data is non- 
vanishing only at the birth points e and / is orthogonal to any evolved pre-birth mode. This 
is because (by the evolution rule) the latter vanish at e and /, and also, by the field equation 
Q6.5| )) the averages on IF(e) and IF(f) vanish. For two of the birth modes we thus take the 
symmetric mode S'i(e) = Si(f) = 1 and the anti-symmetric mode A\(e) = —A\(f) = — 1, 
where the initial data vanishes at all other points. 

For the other two birth modes we set the data to zero at the birth time. The symmetric 
mode with 82(b) = 1 (and zero elsewhere) is orthogonal to all pre-birth- modes since the 
point b is not to the future of any point except e and /. The fourth and final birth 
mode is the anti-symmetric mode defined by ^(a) = — ^(c) = —1, A2Q)) = 0, and the 
requirement that A2 alternates between +1 and —1 on the other points around the slice. 
This is orthogonal to any pre-birth-mode since the sum of the values to the future of any 
birth time point other than the e and / vanishes. 

Being real, the four birth modes Sip and Ai^ have vanishing norm, and the non-zero 
inner products between them are (Si, £2) = (Ai,A 2 ) = i- We seek two complex linear 
combinations that have positive norm and are orthogonal to both of their conjugates. The 
two modes 

S = (St - iS 2 )/V2, A = (A\ — iA 2 )/V2, (6.6) 

have positive-norm, and satisfy the orthogonality requirement, since all modes are orthog- 
onal to their own conjugate, and symmetric modes are orthogonal to antisymmetric ones. 
Other choices than Q6.6| ) could be made. Like for the zero-mode ( |5.14 ), the factor of i 



could be replaced by £7 (and the denominator replaced by \/2Hej). Moreover, since there 
are two independent positive-norm modes, there is further freedom that involves mixing 
between them. The different Fock vacua so obtained are called two-mode squeezed states. 
For the present model we shall make the simple choice given by ( |6.6| ). 

6.4 Growth-induced particle creation 

Now we have all the ingredients for the particle production computation. We wish to 
compute the mean occupation number of each out-mode when the state is the Fock in- 
vacuum. 

The out-number operator is -/V^ out = a \ ou t a fc,°ut- To compute its mean value one needs 
to express the annihilation operator ou t for the out-mode Xk in terms of the in-modes 
£j. We have 

«fe,out = (Xk,4>) (6.7) 
= £(<Xk,6)a< + (Xk,&)4), (6.8) 



- 21 - 



and using the Fock vacuum condition aj|0i n ) = and the commutation relations ( [4.5| ) this 
yields 

<0i„|iV fc ,out|0i„> = ^|(Xfc,e"i>| 2 - (6.9) 



The sum is over the ingoing harmonic modes and zero-mode, as well as the birth modes. 
Formulae for the Bogoliubov coefficients (3ki = (Xki£i) are given in the appendix. 

The spectrum of left-moving cre- 
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Figure 2: Spectrum of created particles on a lattice 
that begins with 999 points and ends with 1001 points 
after a birth process as depicted in Fig. |]. Plotted is 
the Log of the mean number of left-moving particles 
vs. n/10, where the wave number is k = 27rn/1001. 
Black squares give the total, while the other curves 
show the contributions from the left- and right-movers 
(grey and unfilled squares), the symmetric and anti- 
symmetric birth modes (grey and black triangles), 
and the zero- mode (unfilled triangles). The zero- 
mode parameter is 7 = 0.4. 



ated particles is displayed in Fig. g, along 
with the contributions to the mean num- 
ber ( |6.9D from the left-movers, right- 
movers, zero-mode, symmetric birth mode, 
and anti-symmetric birth mode. The 
number of lattice points starts out as 
999, and the zero-mode parameter 7 = 
0.4 in ( |5.14| ) has been chosen. Most of 
the contribution comes from the sym- 
metric birth mode at low wavenumbers, 
from the anti-symmetric birth mode at 
high wave numbers, and from the har- 
monic modes at intermediate wave num- 
bers. The contribution from left-moving 
harmonic modes is several times higher 
than from right-moving ones except at 
very low and very high wave numbers. 
The contribution from the zero-mode is 
far smaller than all the rest. 

The number itself is presumably not 



as relevant as the energy of the created 
particles. As discussed in section 5.2.3, the Hamiltonian — and therefore the energy — is 
not well-defined on a lattice, due to the ambiguity of adding any multiple of 2n/e to 
eigenvalues of the Hamiltonian without altering the evolution operator. Nevertheless, we 
presume that, as a model for the microstructure of spacetime, the degree of excitation rel- 
ative to the translation-invariant Fock vacuum is relevant in determining the gravitational 
back-reaction. We further presume that the smallest positive value of the frequency for a 
given positive-norm mode would determine its gravitating energy, at least for modes with 
frequency much less than 27r/e (mod 2ir/e). 

This energy and its N dependence can be understood roughly from the form of the 
Bogoliubov coefficients in the appendix as follows. The squared Bogoliubov coefficients 
( A.1 , A.2 ) for harmonic in- and out-modes (|5.13[) goes as \Pki\ 2 ~ iV~ 2 (sin k/ sin I). The 
sum over I values gives Xw(V si n = O(NlnN), so the contribution to the number from 
harmonic modes goes as (In N/N) sin k. A given k mode carries an energy k. (Possibly 
we should instead use sinfe, but that difference will not qualitatively affect the result.) 
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Multiplying the number by k and summing over k values yields ^ A: sin A: = 0(N), hence 
the total energy contribution from the harmonic in-mode terms is O(lniV). 

For the symmetric birth mode the Bogoliubov coefficient goes as |/3fcs| 2 ~ (Ns'mk)^ 1 
as long as k is not too close to ir. Multiplying by the energy k then yields the sum 
^2u{k/ sin A;) = O(N), hence the contribution to the energy is 0(1). A similar result holds 
for the antisymmetric birth mode, with a suppression from k near zero. The squared 
zero-mode Bogoliubov coefficient is suppressed by an additional factor of 1/N, yielding an 
energy contribution of order O (In N/N). 

The total created energy due to the birth event is thus of order e~ 1 lnN, where the 
dimensionful constant e setting the scale of length is restored. If every unit of spatial length 
is born accompanied by this amount of energy creation, the result would be incompatible 
with the effective field theory notion that the field remains near the adiabatic ground state. 
A doubling of the size of the universe would produce a Planck energy density of particles 
in this one-dimensional model. 

It is possible as usual to express the out squeezed state in the 

Fock space of the in- modes. If a^i and (3k% are the Bogoliubov coefficients defined by 
afc,out = J2i(. a kiai+Pkia\), this state takes the form exp(-a^, 1 /3fcja]at)|0i n ), where repeated 
indices are summed. In this squeezed state the in- modes are excited and there are specific 
correlations between the excitations of the different in-modes. 

If the ingoing harmonic modes are assumed to be in their ground state, however, no 
amount of fiddling with the birth mode states and the correlations between those and 
the harmonic modes can substantially change the amount of energy created. The out- 
annihilation operator can be expanded as 



Ofe,out — Ofc/ a /,in 



+ Pkl4in + a kBaB + /3fcB%, (6.10) 



where both birth modes are lumped into the one index B. Assuming only that the in-state 
is the vacuum for the harmonic (and zero-) modes, a^inlO) = 0, the mean number evaluates 
to 

(OtalMfc.outlOta) = Yl l^l 2 + W( a k B a B + PkBa ] B )\m 2 ■ (6.11) 
l 



The birth mode state could be chosen so that the second term in ( 6.11 ) vanishes, but this 
would not have a substantial impact on the order of magnitude of the injected energy. The 
spectrum in Fig. || shows that, by themselves, the harmonic modes already contribute an 
order unity fraction. 

Without modifying the assumption that there are no incoming particles, the only 
way to change this conclusion is to modify the evolution rule. This rule affects the in-Fock 
representation, so it affects the state corresponding to the in-Fock vacuum. The expectation 
value of the out-number-operator ( |6.11 ) is sensitive to this via the Bogoliubov coefficients 



Pki = (Xfcj£z)> since the out-modes Xk are n °t defined until the birth time. Some of these 
coefficients can be suppressed by adjusting the evolution rule. One could tune the rule to 
kill individual coefficients, or one could choose a "smoothing" interpolation rule that would 
suppress the contributions from low wave numbers. However, it seems clear that there is 
no way to avoid the 0(e _1 ) energy injection. 
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7. Discussion 

Using the notion of a phase space of half-solutions, we have managed to quantize a linear 
field theory on the background of a lattice with a growing number of points. To ensure 
a well defined symplectic structure on a linear phase space the lattice was assumed to be 
layered by time slices, and to admit at most one backwards evolution, and at least one 
forward evolution from any initial data. We were forced to this rather special structure 
by the requirement of staying very close to standard linear QFT. It may be possible to 
generalize the structure in various ways and still make sense of the quantum theory. One 
obvious small generalization would be to a different type of lattice that is still layered. 

The computation of particle creation in the two-dimensional model showed that the 
addition of two points necessarily injects an energy of order e _1 into the field. If more 
points are added at different times, a similar energy is injected for every unit of length e 
added. Since the in-state is no longer the in-vacuum, the particle production would not 
be identical, but it must clearly be of the same order of magnitude. This would produce 
an enormous energy density of order e -2 during expansion of a model two-dimensional 
universe. Extrapolating to our four-dimensional universe, this behavior is not observed. 
One might hope that in four dimensions the effect of adding points at the lattice scale would 
not be so brutal, but this seems highly unlikely, as the form of the Bogoliubov coefficients 
should be comparable to those in the appendix. 6 This should be expected, since the birth 
process violates time and space translation symmetry on a time and length scale of e, hence 
it it should be accompanied by a violation of energy conservation of order e _1 . 

A way out of this energy catastrophe might arise if many points are added at the same 
time. Then many modes are born at the same time, and wave interference could perhaps 
suppress the net energy injected. This question can be studied in the two-dimensional 
model. 

If interference effects cannot reduce the injected energy, then one would conclude that 
QFT on a background growing lattice is not the right setting to describe mode birth in 
a growing universe. If there is nevertheless a physically correct formulation, it would 
presumably operate at the level where the spacetime itself is treated as dynamical, rather 
than as a fixed background. The complete dynamics might ensure that the system remains 
near the adiabatic vacuum. Perhaps this scenario could be explored using the causal 
dynamical triangulations formulation of quantum gravity [14]. 

Apart from trying to address the fundamental problem of the cosmological vacuum, 
our formulation of mode birth could be applied to field propagation outside a black hole. In 
a geodesic normal coordinate system the spatial metric grows with respect to the geodesic 
proper time, by an amount that depends upon the radius. Discretizing the radial coordinate 
and keeping time continuous one obtains an expanding lattice version of the black hole 
spacetime in which the number of degrees of freedom is constant. The Hawking effect on 
such a falling spatial lattice was studied in It was found that the outgoing modes 

arise via a process analogous to Bloch oscillation, and the continuum Hawking effect is 
recovered when the lattice spacing is small compared to the inverse surface gravity of the 

6 Perhaps the extra IniV would go away in higher dimensions, but that is not our first concern. 
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horizon. Rather than keeping time continuous and allowing the lattice to expand, one could 
instead discretize the time and allow for points to be added so as to keep the lattice scale 
roughly constant. This would bring into play the sort of indeterminate evolution and mode 
birth studied here. In spite of the energy injection by mode birth, the Hawking effect could 
presumably be recovered provided the Hawking occupation number per outgoing mode is 
greater than the number induced by the birth events. 
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A. Bogoliubov coefficients 

In this appendix we present the inner-products between harmonic out-modes and the var- 
ious negative-norm in-modes that appear in the two-dimensional, growing lattice example 
of Section ||. 

For our purposes — calculation of the value of the number operator — we need only the 
norm of the value of the inner products. Consequently, we can freely choose the over-all 
phase of a given mode. For definiteness, here we set the phase of the out-modes to zero at 
point d (see Figure [l] for point labels). For the harmonic, non-zero in-modes, we have set 
the phase to zero at point h. The left-moving mode, when evolved according to the rule 



of Section 6.3.1, has zero phase at point d, the left child of h. The right-moving mode has 
zero phase at point g, the right child of h. We take the zero-mode and the birth modes 
exactly as defined in Section |6.3.2 . 



The formulas correspond to a left-moving harmonic out-mode Xk (we do not consider 
the out-zero- mode) . The results for a right-moving mode are identical up to a phase, except 
that Eqn. ( |A.1| ) now corresponds to the inner product with a right-moving in-mode, and 
Eqn. ( |A.2| ) to that with a left-moving in-mode. All formulas hold for either sign of k; 
i.e. for both positive- and negative-norm out-modes. The positive-norm modes correspond 
to those for which sin(/c) > 0, and vice-versa for the negative-norm modes. 

For the inner-product with a left-moving in-mode of wave- number I, £, , we have: 



N-3 

(tf.Xfc) = 4v_ 2 4v sin + ^ e m <") 



n=0 



, ... M 1 _ e *(JV-2)(l-fc) 
AL At sin H e 1 '" fc + - - 



(A.l) 



with A k N = l/y/N\ smk\. For a right-moving in-mode with wave-number I, we have: 

<£f, Xk) = -e l ^A l N _ 2 A k N sin k. (A.2) 

The above hold for both positive- and negative-norm in-modes, which one can classify by 
the sign of sine of the wave- number. 



- 25 - 



For the negative-norm in-zero-mode, £ , with free-parameter 7, we have: 



A k 

N 



((-1 + ij) sin k + 2 7 e ifc/2 cos §) 



V V2 7 (iV-2) 

For the negative-norm, symmetric birth mode 5, we have: 

(S, Xk) = e~ l2k (2i cos 2 (fc/2) - e ifc / 2 cos(fc/2)) . 

For the negative-norm, anti-symmetric birth mode A, we have: 

(A, X k) = e- i2k (^|) (sin k + ie ifc / 2 sin(fc/2)) . 



(A.3) 



(A.4) 



(A.5) 



To obtain the inner-products with the positive- norm birth or zero- modes, one can use the 
fact that, from the definition of the inner-product, (<f>,ip) = — ((f)* , ip*}* , where * denotes 
complex conjugation. For example, (S,Xk) = — {S,X-k)*- 
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